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Abstract 

We introduce a stochastic dynamics related to the measures that arise 
in harmonic analysis on the infinite-dimensional unitary group. Our dy- 
namics is obtained as a limit of a sequence of natural Markov chains on 
Gelfand-Tsetlin graph. 

We compute finite-dimensional distributions of the limit Markov pro- 
cess, the generator and eigenfunctions of the semigroup related to this 
process. 

The limit process can be identified with Doob h-transform of a family 
of independent diffusions. Space-time correlation functions of the limit 
process have a determinantal form. 

Introduction 

The present paper originated from harmonic analysis on the infinite-dimensional 
unitary group. Decomposition of natural representations of U{oo) into irre- 
ducible ones leads to a family {Pz,w} of probability measures, which depend on 
the two parameters z and w. These measures live on the infinite-dimensional 
domain fi. The definition of the measures can be extended to even lager set of 
values of the parameters. According to this, we use four subscripts z, w, z' , w' , 
instead of two. 

The measures Pz,w and Pz,w,z',w' were introduced in the papers [Olsh] and 
[BOj . Our goal is to construct and study stochastic dynamics related to the 
measures Pz,w,z',w'- 

The structure of measures Pz,w,z',w' substantially depends on whether pa- 
rameters are integers or not. In the present paper we consider the former case, 
the parameters z and w will be integers. Denote these integers by p and q, 
respectively. In our case the support of the measures is a finite-dimensional 
subset of 17 and can be identified with X = [0, The probability dis- 

tributions Pp.q.z',w' were explicitly computed in [B0| . They are given by the 
Jacobi orthogonal polynomial ensemble. 
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We construct a family Jp^q^z',w' (t) of stationary Markov processes in X. Each 
Pp,q,z',w" serves as an invariant distribution of the corresponding Jp^q^z' ,w'(t). 
The processes Jp,q,z' ,w'{t) are obtained as Umits of certain Markov chains on the 
Gelfand-Tsethn graph. These Markov chains arise in a natural way, due to the 
approximation of the infinite-dimensional unitary group by the increasing chain 
of the groups U{N). We call these chains "up-down" chains. Similar Markov 
chains have already appeared earlier, they were studied by Fulman [Ful] . |Fu2] . 
Borodin and Olshanski |B02] and Petrov [Pej 

The proof of the convergence of our Markov chains on the Gelfand-Tsctlin 
graph is based on the special determinantal form of the transition probabilities 
of these Markov chains. We reduce convergence of the transition probabilities 
to convergence of certain matrices. Passing to the limit in these matrices is 
simplified by the fact that we are able to diagonalize them. 

The limit Markov processes Jp,q,z' ,w'it) turn out to have some interesting 
properties. One proves that Jp^q^z' ,w'{t) is a time-dependent determinantal point 
process. It means that its dynamical (space-time) correlation functions have 
determinantal form and can be expressed through the minors of certain extended 
kernel. 

We also explain that Jp^q^z' ,w'{i) can be identified with Doob /i-transform 
of p -I- g -I- 1 independent random motions in [0, 1]. 

We fully describe processes Jp^q^z' ,w'{i)i i-e. we compute their transition 
probabilities and write down generators and eigenfunctions of the Markov semi- 
groups corresponding to the processes. 

The study of the harmonic analysis on the infinite-dimensional unitary group 
shows numerous connections with infinite symmetric group S'(oo) (see jKOVj 
and [BO. Part (m) in Intoduction] ). The constructions of the present paper 
are similar to the constructions of |B02| . where the dynamics related to S(oo) 
were studied. 

It turned out that Markov chains of the present paper have lots of similarities 
with the ones considered in [Gor] . We use some ideas and formulas of [Gor| . 

We want to emphasize that while in our case of integral parameters all 
processes live in a finite-dimensional space, in the case of arbitrary parameters 
the space becomes infinite-dimensional. Thus, our case can be viewed as a 
degeneration of the general case. The problem of constructing dynamics for 
non-integral parameters remains open, it seems like one has to use different 
arguments for that case. 

The paper is organized in the following way. In the first four sections we 
do some preparatory work: In Section 1 we introduce Markov chains related to 
the Gelfand-Tsetlin graph. In Section 2 we present the measures Pz,w,z',w' and 
formulate the problem, solved in the present paper. In Section 3 we study basic 
properties of the Markov chains under consideration. Finally, In Section 4 we 
express transition probabilities of the Markov chains in the determinantal form, 
which is convenient for limit transitions. 

In Section 5 we state and prove the main results of the paper. We prove the 
existence of the limit process and compute its one-dimensional distributions and 
transition probabilities (see Theorem 15. 1|) . Determinantal form of transition 
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probabilities imply the determinantal property of the limit Markov processes 
(see Proposition 15. 5|) . 

In Section 6 we study Markov semigroup related to the constructed process. 
We compute eigenfunctions and eigenvalues of the semigroup (see Theorem l6.ip : 
then we find a simple expression for the generator of the semigroup (see Theorem 

EM- 

In Section 7 we explain the connection between our process and Doob h- 
transform. 

The author is grateful to G. Olshanski for suggesting the problem and nu- 
merous fruitful discussions. 

The author was partially supported by the Moebius Contest Foundation for 
Young Scientists and Leonhard Euler's Fund of Russian Mathematics Support. 



1 General Markov chains on the Gelfand- 
Tsethn graph 

In this section we introduce Markov chains studied in the paper. 

The Gelfand-Tsetlin graph GT (also known as the graph of signatures) is a 
graded graph, whose vertices are so-called signatures. The A^-th level of the 
graph, denoted by GTat, consists of A^-tuples of integers A = (Ai > • • • > Aat) 
which are are called signatures of the order N. We join two signatures A G GTn 
and II S GT^r+i by an edge and write X ^ fi if 

> Ai > ^2 > ■ ■ • > Aat > flN+l- 

We agree that GTo consists of a single element, the empty signature 0. is 
joined by an edge with every signature from GTi. 

By a path in the Gelfand-Tsetlin graph we mean a sequence of vertices 

A(n) ^ A(n 1) ^ < A(m), X^ e GT^. 

Denote by Dim(A) the number of paths going from to A G GT^r. 
For II e GTjv+i and A £ GT^v set 

^i^^ I ^ |Dini(A)/Dim(^), if A ^ ^ 
1 0, otherwise. 

It is clear that 

E ^'Ha|a^) = i. 

The numbers p^{X \ /i) are called cotransition probabilities or "down" tran- 
sition function. 

A sequence {Afjv}Af=o.i,...i where Mjv is a probability measure on GTjy, 
is called a coherent system of distributions provided that for any A'^ > the 
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measures M^r+i and Mjq are consistent with the cotransition probabihties from 
GTtv+i to GTjv, i.e. 

p^{\ I ii)Mn+i{ii) = Mn{\) for any A e GTat. 

Define the support of a coherent system M as the subset 

supp{M) = {A e GT : A/(A) 7^ 0} c GT, 

where M(A) means the measure A/jv(A) of the singleton {A} C GTat. 

Given a coherent system of distributions wc may speak about transition 
probabilities or "up" transition functions. The transition probabihties (/i | A) 
are defined for aU A G supp{M) by 

pHt^ I A) = • P^A I /i), A e GTN,fi e GTn+i,Mn{X) ^ 0. 

Any coherent system of distributions {Afjv} defines a Markov chain H{t), t = 
0, 1, . . . on the state set supp{M). H{t) takes values in Lt = GTj P| supp{M) and 
its one-dimensional distributions are given by Mt- The transition probabilities 
of H{t) are precisely the numbers | A): 

Prob{H(t + 1) = /I I H{t) = A} = p^p I A). 

Note that, while the transition probabilities of H{t) depend on the coherent 
system, the cotransition probabilities Prob{i/(t) = A | H{t + 1) = ^} are 
nothing but the quantities p^ (A | /i) , and they depend solely on the structure of 
the Gelfand-Tsetlin graph. 

We call H{t) the up chain, corresponding to {M^q}. 

We also define a family of stationary Markov chains Tiq{t), = 0, 1, . . . . 
The state space of Tff{t) is Ljv = GT ^ supp{M) . The distribution of TAr(t) 
is given by Afjv. The transition probabilities are given by the composition of 
"up" and "down" transition functions, from Lf^ to L^^i and then back to Lj^. 
Denote by p%{X \ A') the transition probabilities of T/y. We have 

P?^(A|A')= J2 pHAlMy(A^IA'). 

We call T/v(t) the N-th level up-down chain, corresponding to {Mn}- 
In the present paper we study the behavior of the up-down chains, corre- 
sponding to certain coherent systems, as N ^ <x. 

Remarks. Our definitions makes sense not only for the Gelfand-Tsetlin 
graph GT. In their papers Kerov and Vershik introduced transition and co- 
transition probabilities for arbitrary branching graphs (see e.g. |K0V1 Section 
9] and references therein). Thus, we can define Markov chains similar to Hit) 
and T^lt) for an arbitrary branching graph. The up-down chains for the Young 
graph were recently studied by Borodin and Olshanski in [BO 2) . Our work was 
infiuenced by that paper. 
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2 Measures Pz,w,z',w' and construction of the limit 
process 

In the paper [Olshj a 4-parameter family of coherent systems Af^'*"'^ was 
introduced. These measures appear in a natural way in harmonic analysis on 
the infinite-dimensional unitary group. They are given by 
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TT 

i=i " ^» + + N + l + Xi- i)Tiz' -X,+ i)T{w' + N +l + X,-i)' 

where z, w, z' , w' are complex numbers, the quadruple (z, w, z', w') belongs to 
the set of admissible values (see |01shl Definition 7.6]) and Sn{z,w, z' ,w') is a 
normalization constant 

o ( I TT V{z + w + z' + w' + i) 

SN\,z, w,z ,w ) JLl ^ ^ _^ ^ ^^Y{z' + W + i)T{z' + w' + i)T{i) ■ 

Let us assume that z — k and w — I, where k,l G "E, k + l >Q. It was shown 
in [Olshj that (fc, Z, z', w') forms an admissible quadruple of parameters if z' and 
w' are real and z' ~ k > ~1, w' — I > —1. 

Note that for any integer n the shift 

k i-^ k + n, I I — n, z' i—* z' + n, w' ^ w' — n Xi i—t Xi + n 

leaves the probability distributions Mjq invariant. This means that essentially 
the coherent system {M^'''^ } depends on three, not four parameters. From 
now on we assume that fc = p > and I = z > p — \ and w > —\. 

It is easily seen from the definition of M^°'^ that supp{M) consists of 
the signatures A e GTjv such that 

P>Xi>--->Xn>0- 

We always assume below that A satisfies these inequalities. 

Set Xn,p = {0, 1, . . .,7V + p- 1}. Observe that 0<X^-i + N<N+p-l. 
Let us associate to A a collection X{X) of distinct points in Xjv_p as follows 

X{X) = Xn.p \{Xi-1 + N,X2-2 + N,...,Xn- N + N}. 

It is clear that |X(A)| — p and ^(A) defines A uniquely. Let 

X{X) ^ {0 < xi < X2 < ■ ■ ■ < Xp < N + p - 1}. 

Remark. Equivalently A i-^ X{X) can be described in the following way. 
First, we identify A with a Young diagram. Row lengths of this diagram are Aj. 
Then we consider the transposed diagram and its row lengths (in other words. 
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we consider heights of the cohimns of our diagram). We arrive to a collection 
of numbers 

< Ml < < ■ • • < Mp < ^• 

Finally, we set Xi — iii + i — 1. Equivalence of this definition with the one given 
above follows, e.g. from |Macl (1.7)]. 

We work with X{X) instead of A by the following reasons. First, |X(A)| = 
p and, thus, it docs not depend on N. Consequently, we can identify X{X) 
corresponding to different N with elements of one fixed space. Second, transition 
to the next level of the Gelfand-Tsetlinn graph is very simple in the ^(A) 
interpretation. We will provide more details later. 

Denote by P^^ the pushforward of the probability measure M^'°'^ 
under the map A i— > -''^(A). It is clear that Pj\f is a probability measure on 
(Xjv,p)^ and its support is 

supp{Pn) = {{xi, ... ,Xp) e iXN^p)P : xi < X2 < ■ ■ ■ < Xp}. 

For any = 0, 1, 2 . . . we embed the set {Xn.pY into X ~ [0, 1]^ as follows: 

Denote by P"^^ '™ the pushforward of the measure P^^ under the embedding 

TTat. 

Proposition 2.1. As N ^ oo the measures P^^ weakly converge to a mea- 
sure t^p,z' .w' ■ 

The measure ^p^z',w' givefi by its density function Pp^z\w' 
Pp.z' ,w' {x\ , . . . , Xp^ 

Bp.z'w' ■ n - ^jY n (1 - X^Y ~'P, Xi < X2 <, . . . ,< Xp, 

i<j i=l 

0, otherwise, 
where Bp^z'w' is a normalization constant. 

Proof. A slightly different version of this proposition was proved in |B01 Theo- 
rem 11.6]. We will also verify this proposition by straightforward computations 
later (see Proposition 15. 2p . □ 

The aim of the present paper is to introduce a stationary stochastic process, 
which has /x^-^ as an equilibrium measure. We construct this process as a 
limit of the up-down Markov chains on the Gelfand-Tsetlin graph. 

Denote by Xp,z'.w'{t) the image of the up Markov process H{t), correspond- 
ing to the coherent system {M^'°'^ '™ }, under the map A i-^ -^(A). 

Denote by U^^, ^i{t) the image of the up-down Markov process TN{t), cor- 
responding to the coherent system {Mfj'^'^ }, under the map A t-^ -'^(A). 
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Finally, set 

Below we prove the following theorem: 

Theorem 2.2. There exists a limit stationary Markov process Jp.z'.w'{t) on X . 
The finite- dimensional distributions of J^^j ^, {t) converge as N ^ oo to the 
corresponding finite- dimensional distributions of Jp^z',w'{t). 

We will compute the transition probabilities of Jp^z',w'{t) and the generator 

of Jp,z',w'{t)- 

3 Properties of Markov chain Xp^z',w'{t) 

In this section we compute one-dimensional distributions and transition prob- 
abilities of the Markov chain Xp_z',w'it)- In what follows wc omit indices and 
write X{t). 



^1 ,,,/ 



Recall that one-dimensional distributions of X{N) coincide with Pf^^ . 
pp,z ,w ^ pushforward of the measure under the map A ^ X{\), 

and 

M^''"'"''"''(A) = {Sn{z,w,z',w'))-^ ■Uim'iX) 

N 



n 



r(z - A, + i)V{w + N + l + i)T{z' - Aj + i)V{w' + N +l + \-i). 

By Weyl's dimension formula (see e.g. [Zhj ): 

„ , , N TT A," — A," + 7 — i 
Dim A = TT — ^—^ 

l<i<j<N ■' 

= n AeCT.. 

l<i<j<N ■' 

Lemma 3.1. Lei ^ = {0, 1, . . . , A:}, X C A, X = A \ X. Denote 

xex,yex 

v{x)= n (y-^)- 

x<y 

Then we have 

x£X '=1 
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Proof. It is clear that 



^wn n(^-*)l =v{A)-v{x). 



To complete the proof we observe that 

k 



i=l 

and 



Y[ix-i) =x\{k-xy.. 



□ 



As a corollary of Lemma 13.11 and Weyl's dimension formula we have 
Lemma 3.2. Let A e GT^ and X{X) = X ^ {xi < X2 < ■ ■ ■ < Xp), ) then 

Now denote 

, . T{z' + N-x)r{w' +X + 1) 

"'^("^ ^ m+p-xnx + i) ■ 

Below we use the Pochhammer symbol 

, , T(a + k) , . . , , 

r(a) 

Proposition 3.3. For X = {xi < • • • < Xp) we have 

P^/'-'iX) = Z^/^-' ■ V\X)f[wN{x,), 



where '™ is a normalization constant 

N ,.s p 

yp,z',w' _ TT l^yp TT 



{z' + w' + i)p f-J^ T{w' + i + l)T{z' -i + 1)) 
Proof. Straightforward computation using Lemmas 13.11 and 13.21 □ 
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Now suppose that A G GT^, fi e GTjv+i, X = {xi < ■ ■ ■ < Xp) = X{X), 
X' = {x[ < ■ ■ ■ < x'p) = X{^). We write X -< X', provided that A -< ^. Note 
that X ^ X' a and only if for every i, x'^ ~ Xi ot x'^ — Xi + 1. 

Remark. Consider an arbitrary path r = (r(0), r(l), . . . ) in GT. Let 
X{N) = X{TiN)) = (xf < ••• < x^). For every 1 < i < p and N > 
we put a point on the plane with coordinates {N, x'j^) and draw a segment 
connecting {N,xf) with {N + In this way every path t corresponds 

to a collection of p non-intersecting paths on the plane Z^. This interpretation 
shows deep similarities between the objects of the present paper and the objects 
of the paper [Gorj where collections of non-intersecting paths inside a hexagon 
were studied. 

Our next aim is to obtain explicit formulas for the transition probabilities 
of the process Xp,z',w'{t)- 

Proposition 3.4. Suppose X — {xi < ■ ■ ■ < Xp) £ {X^^pY and X' = {x[ < 
■ ■■ < x'p) e {Xn+i^pY, then 

PToh{X{N + 1)=X' \ X{N) = X} 

i.'+M+i), vfx) n (z' + N-x,) n iw' + l + x,), X^X' 

i:x'^—Xi i:x'^—Xi-\-l 

0, otherwise. 

Proof. By the definition Prob{X(A^ + I) = X' \ X{N) ^ X} in transition func- 
tion I A), where X' = X{fi) and X = -'^^(A). Lemma lB^ and Proposition l3.3l 
imply that 

zK'i"' ■ y'ix') n v(x) ■ n ■ ft (« + * - 1)' 

i—1 i—1 i—l 

Z^/'-' ■ V^{X) fi w^{x,) V{X') ■ n ,n(^lp_^,y ■f{{N + ^)\ 

i—1 i—1 ^ * i—1 

Using ll]) we obtain the desired formula. □ 



4 Determinantal form of transition probabilities 

In this section we introduce determinantal formulas for the transition prob- 
abilities of the Markov chains under consideration. These formulas are very 
important for our arguments. 
Denote 

cf = (1 



N I \ w' + z' + N + 1 
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Let VN{x,y) be {N+p) x {N+p-1) matrix (0 < a: < N+p-1, 0<y < N+p) 
given by 



{w'+x+l){x+l) 
(x+N){w'+z'+N+l) ' 



Vn(x,v)^< / iz'+N-x){p+N-x) 
I V (p+2:)(t«'+2'+A^+l)^ 



y ^ X +1 



y = x 



^0, otherwise. 



Proposition 4.1. 

Prob{X(iV + 1) = X' I X(A^) = X} 



' ' dct[wAr(x,,x'0]jj = l,...,p • — ^j— rr, (2) 



'pp/'^^{x) n 

where [vN{xi, x'^)]i_j=i^,,,^p is a submatrix of the matrix VN{x,y). 

Proof. Observe that VN(x,y) is a two-diagonal matrix. Any submatrix of a 
two-diagonal matrix, which has non-zero determinant, is block-diagonal, where 
each block is either upper or lower triangular matrix. Thus, any non-zero minor 
is a product of suitable matrix elements. Consequently, if X X', then 



. ' ■ det[vNix^,x'j)],j=i,...,p ■ 1 

PP/' ''" {X) ^^=0 Cf 



X 



' (;z^ + iV-a;,)(p + iV- xO -pp / (u/ + a; + + 1) 
(p + Xi)(w;' + z' + 7VH 

r(z'-Hw' + iV+ 1) V{X') 



n yz' -t- iV - J 1^;^ -I- iV - J y-r 
V (33 + Xi)(w;' + z' + 7V + l) \ (x + N)(w' + z' + N +1) 



~Tiz'+w'+p+N+i)vix) n (-'+^--^) n 

= Prob{X(A^ +1) = X' \ X{N) = X}, 

and both sides of ([2]) are equal to zero if X 7^ X'. □ 

Let Q^, ff^p_i{x) be the Hahn polynomial of degree k. These polyno- 
mials are orthogonal with respect to the weight wn{x) (see |KSj ). 
Denote 

Qt',z'^pM+p~li^)V'U^Nix). 



\J (Qw',z'-p.,N+p~l^ Qw',z'-p,N+p-l) 



where iQt'.z'-p,N+p-i^Qt'.z'-p,N+p-i) is the squared norm of Qt,^,,_p^N+p-i 
in L2({0, 1, . . . , TV +p — 1}, wn{x)). We provide the exact value of this norm in 
Appendix. 
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Functions f^, k = 0,1, . . . , N + p — 1 form an orthonormal bases in 
L2({0, 1, ■ ■ ■ , N + p ~ 1}) (this L2 is with respect to the uniform measure). 

Remark. Here and below we use definitions of classical orthogonal polyno- 
mials given in [KS| . Thus 

k f-k,-x,k + a + (3 + l\ 



-M, a + 1 



We also use the book [NSUj . which uses slightly different definition of these 
polynomials (two definitions differ by a factor). All formulas in the present 
paper are written according to the definitions of (KSj . 

Proposition 4.2. We have 

N+p-l 
i=0 

Equivalently, 

vn ~ (Pn)'^ ■ Cn ■ (-Fjv+i), 

where vn is {N + p) x {N + p + 1) matrix vn{x, y), F„i is {m + p) x {m + p) 
matrix fmij)i,j=o,...,m+p-i and Cn is {N + p) x (N + p + I) matrix given by 



cf , ifi = j, 
0, otherwise. 



Proof. The argument repeats |Gor| Lemma 8]. More details are given in Ap- 
pendix. □ 



Our next aim is to obtain a determinantal representation for the transition 
probabilities of the Markov chain 11^^, ^,{t). 

Lemma 4.3. If X' < X then 



Prob{X(iV) X' I X{N + l)=X} 

P^''-^'{X~) 

det[vN{x'„Xj)]i,j=i^...,p ■ 1 
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Proof. 

Prob{X(Ar) = X' I X{N +l)=X} 

= Prob{X(7V + 1) = X I X^N + 1) = X'}, P'^^^^^^^) = 



Prob{X(Ar + 1) = X} 



• det[iiAr(x-,a;j)],j=i,...,p 



- , , - • det[wjv(Xi,Xj)]ij=i,...,p • 1 

p^'ir (^) nLo cf 



JV+l 

□ 

Proposition 4.4. Transition probabilities of the process U^.^,^^,{t) are given by 



Prob{[/p^,,,^,(t + 1) = X' I C/p^,,,^,(t) = X} 

, 1 
= =det[Mjv(x.,Xj-)],,j^i,...,p • 1 ^ , 

^pp,^',^'^x) nLo(cf)' 

w/iere [ujv(a;j,a;^)]ij=i,...,p is a submatrix of the {N + p) x {N +p) matrix 

UN = VN ■ (fAf)^ 

Proof. By the definition of U^^^,^^, we liave 

Prob{C/p^,,,,,(t + 1) = X' I C/p^,, = X} 
= Prob{X(Ar) = X' I X{N+l) = y}-Prob{X(Ar+l) = Y \ X{N) = X}. 

Y 

Thus, 

Prob{f/p^,,,^,(t + 1) = X' I U^^,,^^,{t) = X} 

-^■v " '" ("^') del [rrs,: {■i:,,!Jj)],.j^i,...,ij del [r^y (.i;-, <y^)],,j = i....,p 



^p^'^''"''(x) yi<y2<...<yp nr=o(cf) 

Let ^ be p X (A^ + p + 1) matrix given by 

Aij = VN{xi,j - 1), i = l,...,p, j = I, . . . ,N + P+1. 
Let B he {N +P+1) X p matrix given by 

Bij =VN{i-l,Xj), i = l,...,N + p+l, j = l,...,p. 



12 



Denote by A'^^^---^'' the square submatrix of the matrix A consisting of the 
columns ji, j2, • ■ • ,jp] denote by B-Ju^-jp ^j^g square submatrix of the matrix 
B consisting of the rows ji, j2, ■ ■ ■ ijp- By Cauchy-Binet identity 

l<ji<j2<---<3p<N+p+l 

Observe that {AB)ij — UN{xi, x'^). Consequently, 

Prob{[/p^,,,„,(t + l)=X'\ = X} 

__ \l Pn''''''" i^') det[uAr(a;i,x^-)]jj=i,...,p 

□ 

Proposition 4.5. Let k e Z+, then 

Prob{[/p^,,,„,(t + k)^X'\ f/p^,,.„, W = X} 

_ \/ -P^'^''"" (^0 det[wAr,fc(a;i,a;^-)]ij=i,...,p 

^pp'''''^'{x) nr=o(cf)''= 

Proof. We have 

Prob{f/^^,,,„,(t + fc) = X' I = X} 

= J2 Proh{U^^,,^^,{t + k)=X'\ U^^,,^^,{t + 1) = n 

Y 

■ Prob{(7p^,,,„,(t + 1) = r I - X} 

= ^ Prob{{/p^,,^^,(t + fc - 1) = X' I C/p^,,,^,(t) = n 

Y 

■ Prob{[/p^,,^^,(t + 1) = y I ;7,^,,,„,(i) - X}. 

Applying Proposition l4.4l and Cauchy-Binet identity by induction we obtain the 
formula ([3]), where [wN,kixi, x'j)]ij=i^,,,^p is a submatrix of the matrix WN,k = 

(ujv)'^. Consequently, 
Recall that 

{F,nh ^ fraij)^ i , j = , . . . , 771 + p - 1. 
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Since functions form an orthonormal basis, Fm ■ {F„i)'^ = Id, where Id means 
the identity matrix. Thus, 

This equahty imphes Proposition [531 n 

5 Limit process 

Recall that we are interested in studying the stationary process 

Jp,z',w-{t) = lim Jp,,,_yj,{t), 



where 
and 



ttn : {xi,X2, ■ 



Let us introduce some notation. Denote by Wp^z',w'ix) the density function 
of i?-distribution: 

Wp,2',«,'(a;) = a;"''(l - xY'-p, < a; < 1. 

Let Jac'^,_p^,{x) be the orthogonal Jacobi polynomial of the power k. These 
polynomials are defined for x € (0, 1) and are orthogonal with respect to the 
weight function Wp^z' ,w' ■ According to |KSj 

(a + 1)„ (-n, n + a + [3 +1 



Note that while Jacobi polynomials are often defined for x E (—1,1), in the 
present paper we scale and shift the domain of definition, and our polynomials 
live on (0, 1). 
Denote 



jp.z' ,w' (•^) 



where {Jo,c'^'-p w'l Jo^c'^'-p w') ^^'^ squared norm of Jo,Cz'-pw' 
L2([0, l],Wp^z' ,w')- We provide the exact value of this norm in Appendix. 

It is clear that functions jp z' w'j k = 0,1, . . . form an orthonormal system 
in L2([0, 1]) with Lebesgue measure. 

Finally, set 

Kp,z',w'{i) = i{i + w' + z' + 1 - p). 

In this section we prove the existence of the limit process Jp^z'.w' (t) (Theorem 
I2.2p and compute its one-dimensional distributions and transition probabilities. 
Our main result is the following. 
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Theorem 5.1. Finite- dimensional distributions ofJ^^, ^,{t) weakly converge as 
N ^ oo to the corresponding finite- dimensional distributions of a limit process 
Jp,z',w'{t)- Jp,z',tu'{t) is a stationary Markov process. Its initial distribution is 
given by the density 

p 

Pp.z' . w' (^) — I^p.z' . w' {-^i -^j ) ^p^z' . w' (-^i ) ■ 

The transition probabilities are given by the density 



V Pp,z',w'yX) 

where 

oo 

•^p.z' ,w'{-'^tV) ~ ^ ^ ^ '''' ^ \'}p,z',w'i-'^)jp,z',w'iy) 
i=0 

and 

p-i 

^p,z' ,w' — ^ I^p^z' ^w' {}) • 
1=0 

Remark. Determinantal form of transition probabilities implies that 
Jp,z',w'{t) is space-time (dynamical) determinantal process. See Proposition 
15.51 for the exact statement. 

Let us begin the proof. 

In the first place we want to recompute the density of the one-dimensional 
distribution of Jp^z',w'{t) (it was originally computed in [B0| ). 

Proposition 5.2. We have 

PToh{Jp^z',w'{t) e dX} = Pp^z' ,w'{X)dxi . . .dxp 

p 

I^p,z\w' J_('^i '^J ) -^i (1 ^dXij 

i<j i—1 

where X = {xi < X2 < • ■ • < Xp) 

Proof. Clearly, one-dimensional distributions of J^^, ^,{t) are 

PT0h{J^^,^^,it) = X} = P^-'^-'((7V+p- 1) -X), 

where X =^ {xi < X2 < ■ ■ ■ < Xp) and {N + p - 1) ■ X = {{N + p - l)xi < 
{N + p- l)x2 <■ <iN+p- l)xp) e {0,1,...,N +p- 1}P. 
If the limit below exist, then the following equality holds 

Prob{Jj,,,-,^.(t) e dX} = \im {N+p- If ■ Prob{J„~,, „,(i) - X}. (4) 
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Applying the formula 



r(M + b) 

(here A ^ B means lim ^ = 1) we get: 

Prob{ J,^,, ^X}^ Z^' . X\{{N+p^ l)x, ~{N+p- l)x, f 

1^ T{z' + N-{N+p- l)x^)^(w' + {N+p- 1)X^ + 1) 
■ TiN + p-{N + p- l)x^)T{{N + p- l)x, + 1) 

~ Z'/'""' ■ l[{x. -x, f-{N+p- l)^'(f-i) 

i<j 

P 

■ \{{{N +p-l)^{N+p^ l)x,y'-P{{N +p- l)x,r' 
1=1 

= Z^/^^' -{N+p- . _ ^^)2 . -Q(^ _ x,r'-Pxf 

i<j i—1 

Substituting from Proposition 13.31 we obtain 

□ 

Remark. Note that the convergence in ^ is uniform on any compact set 

D C {{xi,X2, . . . , Xp) : < xi < X2 < ■ ■ ■ < Xp < 1} . 
Next we concentrate on the multidimensional distributions of Jp.z'.w'{t)- 
Proposition 5.3. Consider n distinct times 

<ti <t2 < ■ ■ ■ <tn. 

As N tends to infinity, n -dimensional distribution of J^^, ^,{t) corresponding 
to t = ti, . . . ,tn converges. The limit distribution is given by its density 

p{ti,0 < x\ < ■ ■ ■ < xl < 1; . . . ;i„,0 < x^ < ■ ■ ■ < Xp < 1) 



\J Pp,z' ,w' . . . , xjj) 



det[J-;,V;^,(xi,x2)],,,=i,2,...,, det[J^:^7y{xr\x])h,j=i. 



2,...,p 



X \J Pp,z' ,w'ixi , . . . , Xp), (5) 
where J^p^' w'{x,y) and Kp^z',w' are the same as in Theorem \5.1l 
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Proof. Recall that J^^, „,,(t) is a Markov process. Proposition [23] gives its one- 
dimensional distributions while Proposition l4. 51 provides transition probabilities 

Consequently, the joint distribution of J^^, ^i{ti), . . . , J^^, ^/(tn) is given 

by 



Prob{Jp^,, ^,(ti) = X\... = X"} 

det[wN,it,.N^i-lt^.N^]iiN+p-l)x},iN+p~l)x])] 



'?)ji,i=l,...,p 

X 



nr-0 (cf )2 1*2 • ^2 J _2 [ti J 

det[wA,,Lt„.w2j-Lt„-i-JV2j((iV+p- l)x^-\(iV+p- l)a;p],j=i,...,p 



Proposition 15.21 implies that as — > oo, 



^Jpr'''"'{{N +p- 1)X) ■{N + p-l)P^ ^p,.,,.^, {X). 
Next observe that 

\_t-N^\ 

m (cf )2L*-^^J = hm ffl-^]fl + ]] 

N- 



hm (cf )2L*-^'J = hm ((l —\ ( 1 + ]\ 



N~.oo\ {p + N){w' + z' + N +1) ) 

It remains to prove that for t > s, as N ^ oo 

(N+p- lW,Lt.Ar2j_L..iv2j {{N+p- l)x, {N+p- l)y) ^ J^;J^^,{x,y). 

In the following computations we omit parameters and write instead of 
Qw' z'-p N+p-i- ^1^0 write a; instead of {N + p — l)x and y instead of 
{N + p— l)y. We have 

N+p-l 

4=0 

^wp.,M-)^p,,My){N+p-ir'^^'-^- (cf)^L*-A'^j^2Ls.A.^jQ|i)|W 



1=0 



^Wp,,>^^\x)wj,^,,,^.{y){N +p- ir'+''-P{A + B), 
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where 



and 



Now fix large enough / and send N to oo. First, let us examine A. We have 
(see [IHl Section 2.5]) 



Qlo',z'--p,N+p-liiN +P-1)X) 



and the convergence is uniform in x belonging to any compact subset of (0, 1). 
Here we used the relation 

Jacj.,^(l - a;) = (-l)Va4^„(a;). 

By straightforward computation one proves that 



lim 



N^oo iQl„,^z'-p,N+p-l^Qw'.z'-p,N+p-l) iJ^(^l'-p,w'(^)^ J^'^\'-p,w'i^)) 



in slightly different form this claim was also proved in |NSU| . 
Consequently, 



lim {N + p-l)- Jwp.z>.^,{x)wp,,,.^,{y){N +p- l)^ A 

N—>-OG 



oo 
i=l+l 

Observe that for an arbitrary e > we have 

oo 

if I is large enough. Indeed, 



and 



{x)i 
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where c{x) is a bounded function for x belonging to any compact subset of 
(0, 1). This estimate fohows, for instance, from the recurrence relations on 
Jacobi polynomials. 

Thus, to finish the proof we should show that 



limsup Jwp,.',^>ix)wp,,,,^,iy){N +p- 1)^'+''+'-p\B\ < £(/), 

and e{l) tends to zero when I tends to oo. 
We need the following lemma: 

Lemma 5.4. For large enough N and any i < N + p we have 

where constant c does not depend on N . 
Proof. 

/^^^2^-| iji + w' + z' -p+1) f_ 
^""^ ^ {p + N){w' + z' + N+l) "^N^' 

If I < j^^f^ and N is large enough, then 



c^i In" i 



If i > then 



In(iV) \ in(jv) 

HN)J j 

□ 

Recurrence relations on Hahn polynomials (see e.g. |KS| (1.5.3)] imply that 

\Qw' ,z'-p,N+p-li^)\ < 

where constant c' does not depend on either iV or x. 
Next note that 



(N + p-l)'^'+^'+^-P 

iQw' ,z'-p,N+p-l' Qw' ,z'-p,N+p-l) 
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is bounded from above by e'^ ■«in(«) 

For large enough N and I we obtain the following estimate 

{N +p-ir'+''+^-P\B\ ^ 

N+p-l 
i=l+l 

{N + p- l)«''+2'+i-p 



(O^ , , 

w' ^z' —p. 



AT+p-l' ^u;',z'-p,A'+p-l - 



^ g-ciln^(i)+2c'iln(i)+c"iln(i) ^ g-ciln^(i) ^ ^-ci _ ^ 

i=;+l i=;+l i=l ^ ^ 

— cL 

Observing that ^_si when ^ ^ oo, completes the proof of Proposition 

□ 

Since multidimensional distributions of the process J^^, (t) converge, we 
define Jp^z'.w'{t) as a limit process. Formulas ([S]) for the multidimensional dis- 
tributions imply that Jp^z'.w'{t) is a stationary Markov process with initial dis- 
tribution and transition probabilities as in Theorem 15.11 Thus, Theorems 12.21 
and 15.11 are proved. 

The Markov process Jp,z'.w'{t) has one interesting feature. Its dynamical 
(space-time) correlation functions can be expressed as minors of a certain ex- 
tended kernel. 

Let Pn{xi, ti;x2,t2;--.;Xn,tn) be the nth correlation function of Jp^z',w'{t)- 
Informally p„ can be defined by 

Prob{xi €: Jp^z' ,w' {pl) 1 • ■ • 1 "^n ^ ^p,2',iij'(^n)} (-^1 ; ^1 : • ■ ■ ) -^n : tji^dx\ . . . dXji ■ 

Proposition 5.5. Consider n distinct points (xi, ti), . . . , (a;„, i„). We have 
w/iere 



Kerp,2/,t„/(a;,<; ?/, s) 



i=0 



We are not going to give the proof of this claim here. It follows from Theorem 
Oand Eynard-Metha theorem. See 'EM^ and ^603} Section 7.4]. 

^ is a good survey on determinantal point processes. Additional informa- 
tion about extended kernels can be found in e.g. [EM| . [j] . |TW| . 
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6 Markov semigroup of the limit process 

Let Wp be the Weyl chamber. 

Wp = {{xi,...,xp) e [0,l]f :xi <X2 <---<Xp}. 

Clearly, Wp is a state space of the process Jp.z' ,w'{t)- Recall that fip^z',w' is the 
one-dimensional distribution of Jp^z' ,w'{t)- According to Theorem 15.11 ^ip^z'.w' 
is a probability measure on Wp given by its density Pp^z'.w'- 

Denote by M^ .,, ,^, the Markov semigroup of the process Jp,z',w'{t) on 
L2{Wp,iJ,p,z',w')- M^p^^'^w' acts on a function f{X) £ L2{Wp, fip^z' ,w') as fol- 
lows 

(M^,,,.^,/)(X) = / F{Y)Pl,,^^,{Y I X)dY, 

where the integration is performed with respect to the Lebesgue measure on 
Wp. 

In this section we are going to compute the eigenfunctions and the generator 

Recall that Jac^,_p (x) is Jacobi polynomial of degree k on (0, 1). Now we 
want to introduce multi- dimensional Jacobi polynomials. Let A — (Ai > A2 > 
• • • > Ap > 0) be a partition of n. Denote 

. A , , _ det[Ja4'+P~!(xj-)],^j^i,...,p 

^^,{xi, . . . ,Xp} — Y{{xi-x ) 

i>j 

It is clear that Jac^,_^ ,^,{xi, . . . ,Xp) is a symmetric polynomial in p variables 
of degree n = Ai -I- • • • + Ap. 

We may view ^ac^,_p_^,(xi, . . . function on Wp. In what follows 

we use the same notation for these functions. 

It is well-known that the polynomials Jac^,_p ^i{xi, . . . , Xp) form an orthog- 
onal basis in L2{Wp, fip^z',w')- 

Theorem 6.1. Polynomials Jac^,_p ^,{xi, . . . ,Xp) are eigenfunctions of the 
Markov semigroup M* ^, ^, . 

Ml,z',n,'{Jacl,_p^^,) = c{\,t)Jac^^^,^^,, 

where 

,^ , £ ^p,.',™'(p-*)--ffp,.',™'(A.+p-*)') 

c{\,t) = e ^'=1 ^ 
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Proof. 



i.Ml,z',w'J0'C^z'-p,w'){xi,X2, ■■■,Xp) 



A , N Vpp,^', wivi, ■■■,yp) 



'yi< - <yp y Pp,z' ,w' (xi, ■ ■ ■ , Xp) 

X e*^p.-'.™' • det[J^.^,^^,{xi,yj)]ij=i^2,...,pdyi ...dyp 

■/ Jac^,_p,^,{yl,...,yp)■\[{y^ 



nl^ \ I """^z'-p,w\yi-T ■ ■ ^yp/ livyi Vj) 



i>3 

Y{i=i,pVwp,^.^^,{yi) 



^'"^ -<:lc^[Jp^z',w'{xi^Vj)]i,j=i,2,...,pdyi---dyp 



Tii^l^p V^p,z',w'i^i) 

Note that the integrand is a symmetric function in yi,...,yp. Thus, 
we can integrate over p-dimensional cube (0,1)^, instead of the simplex 
< yi < • • • < j/p < 1. We arrive to the foUowing expression 



det[Jac^;+^ 



xdet 



dyi... dyp. 

-I ij = l,...,p 



Multiplying two matrices under determinants and using orthogonality relations 
for the Jacobi polynomials we obtain the desired formula. □ 

Let Gp^z',w' be the infinitesimal generator of the Markov semigroup ^, 
Theorem 16. II implies the following proposition. 

Proposition 6.2. Polynomials Jac^^^, ^,{xi, . . . ,Xp) are eigenf unctions of 

dp^z' .w' {J<^C^/ _p ,^^, ) — c{\)JaCzr ^p .fjji , 

where 

p 

c(A) = '^Kp^z'.w'ip - «) - ^p,z',tu'(Aj + p-i) 

Our next aim is to obtain explicit formulas for the generators Gp,z\w'- 
Observe that the numbers —Kp^z'.w'{i) ~ —i(i + w' + z' + 1 — p) are eigen- 
values of the Jacobi differential operator corresponding to the eigenvectors 
Jac'^,_p ,^,{x). Namely, denote 

^f'-p,^' ^xa-x)^ + iw' + l-iw' + z' -p + 2)x)-^. 



The following proposition holds (see e.g. [KS[ (1.8.5)]) 
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Proposition 6.3. We have 
Set 

P q2 q 

1^z',w',p = - Xi)-^ + {w' + 1- {w' + z' -p + 2)Xi) — . 

Recall that if X = (a;i, . . . , Xp) then V{X) = H (^i ~ ^j)- 

i>j 

Theorem 6.4. We have 

and the domain of definition of Gp^z',w' consists of all functions 
f € L2{Wp, fip^z' ,w') such that 

'Dz'^w'AVix) ■ f) + Kp^z'^yj'f e L2(Wp, llp,z',w')- 



V(X) 



Remark. After some simplifications the generator Gp^z',w' can be alterna- 
tively expressed as 



Gp,z',w'{f) = 1^z',w',p{f) + X] ~ ^i) H Z 



dXi 



Proof. It is sufficient to verify that 



Gp,z' ,w'iJcic^_pi ^^,) — Y(^x) ' '^^'''"''P (^("^-^ ' '^'^^z-p' ,w') ~^ ■^p,z' ,w' Jac^_p, . 
We have 

1 ~ ' A All.- . 7„„A 



ViX) 



T^z',w',p<iet[Jac^,_^.^,{xj)]ij=i^...^p ^ 

r J^n.z' .w' "J (^^r. — n'm' • 



Let us write the determinant as the alternating sum of products 

A.^\{jac'j^^J,{x,^,^). 

(Here tr e is an arbitrary permutation.) 
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Proposition 16.21 implies that 



Consequently, 



T^z',w',pdet[Jac^,_p.^,{xj)]ij=i^,,,^p ^ 

.(.„.-(t,,,„.,..-,,))«_,,,..,,.,,..,,, 



□ 



7 Doob h-transform interpretation 

In this section we explain that the process Jp^z',w'{t) can be viewed as p in- 
dependent identically distributed diffusions conditioned (in the sense of Doob) 
never to collide, i.e. it is the Doob h-transform of p independent processes. 

Proposition 7.1. We have 

Proof. This claim follows from Proposition 16.21 and Theorem 16.41 but we also 
present an independent proof. 

Lemma 7.2. For arbitrary a,b,c E M. set 

(G/)(X)^f](x2 + ax + 6)^+ f-|(p-2K + c^ ^ 
Then GV{X) = 0. 



. , dxj V 3 / dxi 



Proof. See |K01 Lemma 3.1]. Note that originally there was a misprint in [KOj . 
the minus sign was missing. □ 

Lemma 17.21 implies that 
Vz',^',pViX) = (^-^{p -2)-{w' + z' -p + 2)^ 

^(^-'^^ ]ip-2) + w' + z')v{X). 



2 V 3 
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To finish the proof let us compute explicitly Kp^z' ,w''- 



p—i p— 1 



p.z' ,w' ^ 




P{P - 1) / 2p-l 
2 I 3 




□ 



Now consider p independent diffusions with generators D'^,"^^^^,. We call 
these diffusions Jacobi processes; they can be viewed as solutions of certain 
stochastic differential equations. 

The infinitesimal generator of these p diffusions is precisely "Dz' ^w' .p- Propo- 
sition [TIT] implies that we can construct new Markov process, which is a variant 
of Doob /i-transform of p diffusions. Here function h\sV{X). The infinitesimal 
generator of Doob h transform with h = V{X) is precisely Gp,z',w'- The Doob 
/i-transforms and its properties were studied in numerous papers, see e.g. [D] . 



Thus, our process Jp,z',w'{i) can be interpreted as Doob /i-transform of p 
independent diffusions. In many examples Doob /i-transform coincides with the 
conditional process, given that there is no collision of the components. In our 
case, clearly, there are no collisions, trajectories of Jp^z',w'{t) are collections of 
p nonintersecting paths. However the question whether Jp^z'.w'{t) is precisely p 
independent diffusions conditioned on never having a collision between any two 
of its components, remains open. 

8 Appendix. Some formulas for Hahn and Ja- 
cobi polynomials 

The following formulas can be found in e.g. [KSj . 
Lemma 8.1. For a > —1, (3 > —1, 



[Kcm) . 



M 



r{a + x + l)r(/3 + N - x + 1) 
r(a; + l)r(iV-.T + l) 




' Qa,P,M) — / , 



Qa,l3,M i^)Qa,f},M {^) 



(-l)fe(fc + a + 0+ 1)m+i(/3 + l)fcA:!r(a + l)r(/? + 1) 
{2k + a + (3 + l)(a + l)fe(-M)fcM! 



Ski ■ 
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Lemma 8.2. For x, ?/ G {0, 1, . . . , M), 

^ (-l)'^(A; + a + /3 + 1)m+i(p + IjfeK! 



(a + l),(/3+l)M- 



a;!(M-x)! 

The following Lemma was proved in |Gori Lemma 8] 
Lemma 8.3. 



xQk - 1; a, /?, A/ - 1) + {M ~ x)Qk [x; a, (3, M - 1) = MQk{x] a, (3, M). 
Proof of Proposition Recall that our goal is to prove the following equality 



N+p-l 



V (x+N)(w'+z'+N+i) ' y — ^ ^ 



1=0 



(p+x)(w'+z'+N+l) ' 

0, otherwise. 



First, we express left-hand side of ^ through the polynomials ^ j^j. By the 
definition 

Qt',z'-p,N+p-lix) V^'n{x) 



\J (Qw',z'-p.,N+p-l^ Qw',z'-p,N+p-l) 



Lemma O provides the value of iQtj'^,j^p^N+p-i^Qt',z'-p,N+p-i)- ^^ext, 
we express Qt>^z>^p,N+p through Qt' ,z'-p,N+p-i rising Lemma ESI Applying 



Lemma [8T2] completes the proof. □ 

The following lemma in slightly different form can be found in jKSj . 
Lemma 8.4. For a > —1, (3 > ~1, 

1(1 xrx^Jac^ (x)Jac' (x) - r(fc + a + l)r(fc + /3 + 1) 
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